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ABSTRACT 

ZZ Ceti stars exhibit small amplitude photometric pulsations in multiple gravity modes. As the stars 
cool their dominant modes shift to longer periods. We demonstrate that parametric instability limits 
overstable modes to amplitudes similar to those observed. In particular, it reproduces the trend that 
longer period modes have larger amplitudes. 

Parametric instability is a form of resonant 3-mode coupling. It involves the destabilization of a pair of 
stable daughter modes by an overstable parent mode. The 3-modes must satisfy exact angular selection 
rules and approximate frequency resonance. The lowest instability threshold for each parent mode is 
provided by the daughter pair that minimizes (Sto 2 + 7 2 )/^ 2 , where k is the nonlinear coupling constant, 
8u> is the frequency mismatch, and "fd is the energy damping rate of the daughter modes. Parametric 
instability leads to a steady state if \Su>\ > jd, and to limit cycles if \Su>\ < -fd- The former behavior 
characterizes low radial order (n < 3) parent modes, and the latter those of higher n. In either case, the 
overstable mode's amplitude is maintained at close to the instability threshold value. 

Although parametric instability defines an upper envelope for the amplitudes of overstable modes in ZZ 
Ceti stars, other nonlinear mechanisms are required to account for the irregular distribution of amplitudes 
of similar modes and the non-detection of modes with periods longer than 1,200 s. Resonant 3-modc 
interactions involving more than one excited mode may account for the former. Our leading candidate 
for the latter is Kelvin-Hclmholtz instability of the mode-driven shear layer below the convection zone. 

Subject headings: instabilities — white dwarfs — stars:oscillations 



1. INTRODUCTION 

Within an instability strip of width AT e g ~ 10 3 K cen- 
tered at T c ff sa 1.2 x 10 3 K, hydrogen white dwarfs exhibit 
multiple excited gravity modes with 10 2 <J P <J 10 3 s. 
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Conyective driving, originally proposed by Brickhill (1990, 
), is the overstability mechanism (Goldreich & Wu 
, hereafter Paper I) . Individual modes maintain small 
amplitudes; typical fractional flux variations range from a 
few mma to a few tens of mma. 4 

The nonlinear mechanism responsible for saturating 
mode amplitudes has not previously been identified. We 
demonstrate that parametric resonance between an over- 
stable parent g-mode and a pair of lower frequency damped 
daughter g-modes sets an upper envelope to the parent 
modes' amplitudes. 5 Moreover, the envelope we calculate 
reproduces the broad trends found from observational de- 
terminations of mode amplitudes in ZZ Ceti stars. Our 
investigation follow s pioneering work by Dziembowski & 
Krolikowska (1985) on overstable acoustic modes in 5- 
Scuti stars. They showed that parametric resonance with 
damped daughter g-modes saturates the growth of the 
overstable p-modes at approximately their observed am- 
plitudes. 

This paper is comprised of the following parts. In 
we introduce parametric instability for a pair of damped 
daughter modes resonantly coupled to an overstable parent 



mode. We evaluate the parent mode's threshold amplitude 
and describe the evolution of the instability to finite am- 
plitude. §3 is devoted to the choice of optimal daughter 
pairs. We discuss relevant properties of 3-mode coupling 
coefficients, and the constraints imposed by frequency res- 
onance relations and angular selection rules. Evaluation 
of the upper envelope for parent mode amplitudes set by 
parametric resonance is the subject of §4. Numerical re- 
sults arc interpreted in terms of analytic scaling relations 
and compared to observations. §5 contains a discussion of 
a variety of issues leftover from this investigation. Detailed 
derivations are relegated to a series of Appendices. 

The stellar mode ls use d in this investigation were pro- 
vided by Bradley (1996). Their essential characteristics 
are M* = O.6M , log(<?/ cms -2 ) = 8.0, hydrogen layer 
mass 1.5 x 10~ 4 Af„, and helium layer mass 1.5 x 10~ 2 M*. 

2. PARAMETRIC INSTABILITY 

In this sectio n we i ntroduce parametric instability (Lan- 
dau & Lifshitz 1976 ). We present the threshold criterion 
for the instability, and discuss relevant aspects of the sub- 
sequent evolution. Depending upon the parameters, the 
modes either approach a stable steady-state or develop 
limit cycles. We describe the energies attained by the par- 
ent and daughter modes in either case. Many of the re- 



sults in this section were obtained earlier by Dziembowski 

Astronomy Unit, School of Mathematical Sciences, Queen Mary and Westfield College, London, UK 

2 Canadian Institute of Theoretical Astrophysics, University of Toronto, 60 St. George Street, Toronto, Ontario M5S 3H8, Canada; 
wu@cita.utoronto.ca 

3 Theoretical Astrophysics, California Institute of Technology 130-33, Pasadena, CA 91125, USA; pmg@gps.caltech.edu 
4 1 mma of light variation is approximately 0.1% fractional change in flux. 
5 The g-mode dispersion relation allows plenty of good resonances. 

1 



2 



Wu & Goldreich 



(1982) 



2.1. Instability Threshold 

Parametric instability in the context of our investiga- 
tion is a special form of resonant 3-mode coupling. It 
refers to the destabilization of a pair of damped daughter 
modes by an overstable parent mode. The frequencies of 
the three modes satisfy the approximate resonance con- 
dition lu p w u> dl + , where the subscripts p and d\ , di 
denote parent and two daughter modes, respectively. 

Equations governing the temporal evolution of mode 
amplitudes are mos t conv eniently derived fro m an action 
principle (Newcomb [l962] Kumar & Goldreich |l989| ) . The 
amplitude equations take the form 



dA p 
~dT 

dA dl 

dt 
dA d2 

dt 
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~A P - ibipAp + i— =u p KA dl A d2 , (1) 

~M - iuj dl A dl + i-j=uj dl KA p A* d2 , (2) 

~A d2 -iuj d2 A d2 +i-^=cj d2 KApA dl . (3) 



Here, Aj is the complex amplitude of mode j; it is re- 
lated to the mode energy Ej by \Aj\ 2 = Ej. The jj (> 0) 
denote linear energy growth and dampi ng rates, and k is 
the nonlinear coupling constant (cf. §3.1). Our ampli- 
tude equations diffe r only in notation from those given by 
Dziembowski fll982[ ). 

The instability threshold follows from a straightforward 
linear stability analysis applied to equations (g) and (||). 
The effects of nonlinear interactions on the parent mode 
are ignored as is appropriate for infinitesimal daughter 
mode amplitudes. The c ritical parent mode a mplitude sat- 
isfies (Vandakurov 1979, Dziembowski 1982). 



l8K 2 U! dl U d2 
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25uj 



(4) 



where Su> = ui dl + uj d2 — u> p . 

2.2. Dynamics 
The amplitude equations (0)-(||) have a unique equilib- 



rium solution given by 
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7di + ld 2 ~ 7p , 



2Slu 



7di + ld 2 - 7p 



(5) 
(6) 
(7) 
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<f> = 



'pi 



where the complex amplitude 
0j . Note that equa- 



Here, 

Aj may be written as Aj 

tion (0) is almost identical to the threshold criterion (eq. 

6 The angular dependence is described by a spherical harmonic Yj r 



Q) for 7 p <C 7rf which is the case that concerns us. Here 
7d = (7di + ld 2 ) /2 is the characteristic damping rate for 
the daughter modes. It is also worth mentioning that, in 
this limit, the parent mode energy is independent of 7 p . 

The equilibrium state is a stable attractor for mode 
triple ts wi th \6u>\ > j d , a nd un stable otherwise (Wersinger 



et al. 1980, Dziembowski 



1982| ). Figures |a & 0b illustrate 



these two types of behaviors. Triplets with unstable equi- 
libria undergo a variety of limit cycles. These share a num- 
ber of common features. The parent mode's amplitude re- 
mains close to its equilibrium (threshold) value, with slow 
rises on time scale 7" 1 followed by precipitous drops on 

time scale 7J 1 - The daughter modes' amplitudes stay far 
below their equilibrium values for most of the cycle, but 
peak with amplitudes comparable to that of the parent 
mode for a brief interval of length j^ 1 shortly after the 
parent mode amplitude reaches its maximum value. Dur- 
ing this brief interval, the energy which the parent mode 
has slowly accumulated is transferred to and dissipated by 
the daughter modes. So we see that, independent of its 
stability, the equilibrium state defines the parent mode's 
amplitude. This is confirmed by Figure ^i. We employ 
this result in where we predict upper limits for g-mode 
amplitudes in pulsating white dwarfs. The time-averaged 
energies of the daughter modes are also found to be con- 
sistent with equations (Fig. IJd). 

3. CHOOSING THE BEST DAUGHTER PAIRS 

The discussion in the previous section shows that an 
overstable parent mode's amplitude saturates at a value 
close to the threshold for parametric instability. Although 
each overstable parent mode has many potential daughter 
pairs, the most important pair is the one with the lowest 
instability threshold. This section is devoted to identify- 
ing these optimal daughter pairs, a task which separates 
into two independent parts, maximization of k 2 and min- 
imization of (Sui 2 + 7J). 

3.1. Three-Mode Coupling Coefficients 

The 3-mode coupling coefficient characterizes the lowest 
order nonlinear interactions among stellar modes. A com- 
pact form suitable for adiabatic modes under the Cowl ing 
approximation is derived in Kumar & Goldreich (198E); 



-3(r 1 -i)(v.€)r^ + 2^4^}, 



(9) 



where p is the unperturbed pressure, Ti is the adiabatic 
index, £ is the Lagrangian displacement, the symbol ';' de- 
notes covariant derivative, and the integration is over the 
volume of the star. This expression for k is symmetric with 
respect to the three modes. Note that the displacements 
enter only through components of their gradients. The 
present form is not suitable for accurate numerical c om- 
puta tion. We derive a more appropriate version in §A.l 
(eq. ^l|]). 

Each eigenmode of a spherical star is characterized by 
three eigenvalues n, t, m; n is the number of radial nodes 
in the radial displacement eigenfunction, £ is the spheri- 
cal degree, and m is the azimuthal number. 6 Integration 
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Fig. 1. — Parametric instability studied by numerical integration of equations (^)-(^|)- Solid lines represent parent mode amplitudes 
and dashed lines those of daughter modes. We multiply these amplitudes with re to make them dimensionless. For the top panel 
Slu = uj c i 1 + u)ii 2 — ujp = 2 X 10~ 5 s -1 3> Jdi = 7d 2 = 10 — 6 s — 1 , and the system settles into a steady state. The middle panel shows a 
case with 5u> = 10~ 6 s -1 < £i^)jd 1 = 7d 2 = 10 -5 s _1 , for which the mode amplitudes undergo limit cycles. These panels illustrate the two types 
of behavior discussed in § [2.2| The bottom panel displays photospheric flux variations associated with the case shown in the middle panel. 
The daughter mode energies episodically approach that of the parent mode. Around these times their fractional flux amplitudes exceed that 
of the parent mode because they have smaller mode masses. For both simulations, we take mode periods to be P p = 500 s -1 , Pd 2 /Pdi = 0-9, 
and 7 P = 10 -7 s — 1 . 



over solid angle enforces the following selection rules on 
triplets with non- vanishing k: \£d 2 ~ ^di \ < In < I rfi + ^d 2 , 
tp + i-dx + £d 2 even, and m p = +rrid 2 (see §A.l). These 
selection rules guarantee the conservation of angular mo- 
mentum during nonlinear interactions. The magnitude of 
k is largest when the eigenfunctions of the daughter modes 
are radially similar in the upper evanescent zone of the par- 
ent mode. Radial similarity requires near equality of the 
vertical components of WKB wavevectors, k z , where for 
gravity modes k\ w (N 2 /uj 2 — l)A 2 /r 2 (eq. [A2] of Paper 
I) with N 2 being the Brunt- Vaisala frequency, r the ra- 
dius, and A 2 = £(£ + 1). We take N ~ u> v as the major 
contribution to the peak value of | n | comes from the region 
just above z UJp , the upper boundary of the parent mode's 
propagating cavity. The g-mode dispersion relation appli- 
cable for u> -c 10~ 2 s _1 is u> oc t/n (see Fig. 4 of Paper I). 
We find that k is largest when 

nd 2 Ad 2 ^d 2 



"1 - ^4 



1/2 



(10) 



Taking l p = 1, £d 2 = ld x + 1 ; and evaluating the above 
relation for id x < 10, we locate the peak of n to be at 



rid 1 - n d2 



-0.7n r 



(11) 



Since a fraction n" 1 of the nodes of each daughter mode 
lie in the region above z UJp , the peak has a width of order 
n p when measured in — rid 2 ■ Figure || illustrates the 

7 The factor n p arises because the peak in re has a width [n^ — 



behavior of k for a variety of combinations of parent and 
daughter modes. 



As we show in § A. 2 , the maximum value of |k| is of order 



L 



1/2' 



(12) 



Here L is the stellar luminosity and r Wp is the thermal 
timescale at z UJp . This is compared with numerical results 
in Figure ^. Note that the maximum value of k depends 
entirely upon the properties of the parent mode and not 
at all upon those of its daughters. 

3.2. Frequency Mismatch and Damping Rates 

Because the maximum value of is independent of 
td t and £d 2 , we choose these parameters to minimize 
i^ 2 + l 2 d)- 

Consider an £ p , m p parent mode. For each choice of 
£d 1 , £d 2 j there are of order 1& X n 2 daughter pairs for which 
k is close to its maximum value. The factor £d t arises 
from the freedom in choosing , while the factor n 2 
comes from the width of maximum \n\ at each ld x . Re- 
laxing the value of £d 2 subject to the constraint of the 
angular selection rules increases the number of pairs by 
a factor of order £ p . Now replace £d 1 by a running vari- 
able idi- The number of pairs with id/ < £d x is of order 
£p£'d i n 2 ,- The distribution of the 6u> values of these pairs 
is uniform between and n v LOd 1 /nd 1 .' 7 Statistically, the 
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Fig. 2. — Numerical simulations confirm theoretical predictions for average mode energies. We numerically integrate equations for a 

variety of three-mode systems with different frequencies and energy growth/decay rates. Each point on the figure represents one such system. 
For simplicity, we restrict ourselves to systems for which 7^ = ■y^ = 7^. Time-averaged energies of parent and daughter modes are plotted 
against theoretically predicted values (eqs. |5|-H]) m th e upper and lower panels, respectively. Energies multiplied by k 2 are dimensionless. 
The subscript 'd' in the lower panel is taken to De d\ or 6,2- 




Fig. 3. — Coupling coefficient k in erg -1 / 2 for a DA white dwarf model with T c g = 12, 800 K. From top to bottom panels, we consider 
four i v = 1 parent modes with increasing n p . Coupling coefficients between each parent and its daughter pairs are plotted versus — n^,. 
Solid lines denote daughter pairs with = 1, and others those with = 2, 3, 4, 5, and 9. Mode 3 is chosen to minimize \Sui\ subject to the 
constraint = t v + td x ■ t is seen to depend on the radial similarity of the two daughter modes (hcl. jUI) but not on their spherical degrees. 
The FWHM of the peak in |k| is of order n p . Maximum \k\ occurs for — n^ 2 ~ — 0.7n p (eq. [|f 1|). 

minimum frequency mismatch u> oc £/n, is assumed in going from the first to the sec- 

ond relation for Aw. Our estimate for Aw assumes that 
A ~ UJdl - ~ UJp (13) rotation lifts m degeneracy. If it does not, the minimum 

n di ^p^d 1 n p ^di n p ' frequency mismatch is increased by a factor of ld x - 

Next we describe how 7 varies with u> and £. There 
where the low frequency limit of the dispersion relation, are two regimes of relevance to this investigation. In the 
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Fig. 4. — Coupling coefficients as a function of parent mode's radial order. We consider t v = 1 parent modes, and include all daughter pairs 
that have = 1, ...,9 and = dp + ^d 1 - F° r each mode 2, mode 3 is chosen to minimize |<5a;|. The solid line corresponds to the theoretical 
estimate for the maximum at each n v (eq. flld]). L is the stellar luminosity, and 7" w is the thermal time at z u]p . These calculations are 
based on the same stellar model used for Figuretl 




Fig. 5. — Nonadiabatic energy driving and damping rates as a function of mode period. See Paper II for computational details. We consider 
three stellar models with decreasing effective temperature; r t ^ is the thermal time scale at the bottom of the convective zone. Overstable 
modes are marked by solid triangles. Lines connect points associated with_modes of the same t. For each £, 7 exhibits an initial steep rise 
with increasing mode period (eq. followed by a gradual decline (eq. ]15|). Maximum 7 occurs at the transition between quasi-adiabatic 

and strongly nonadiabatic damping. 



quasiadiabatic limit (cf. §4.4 of paper I), 8 



In the strongly nonadiabatic limit (see §B.l), 



7 1 



— In — oc ui 
■nn 7c 



0.75 £0.2 



(15) 



(14) 



where 1Z is the amplitude reflection coefficient at the top 
of the mode's cavity. The transition between the quasia- 



3 Overstable modes are quasiadiabatic so this estimate for 7 applies to them as well as to damped modes. 
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diabatic and strongly nonadiabatic limits is marked by a 
significant reduction of 1Z by radiative diffusion. The be- 
havior of 7 as a function of u> and £ is illustrated in Figure 
Radial similarity of the daughter modes to which a 
given parent mode couples most strongly (eq.|l(|) implies 

Tdi ~ ld 2 ^~ Id- 

The minimum of (Suj 2 + 7J) is attained for daughter 
pairs that satisfy Suj 2 <J 7J. With increasing n p , the value 
of at which this occurs decreases from a few to unity 
(see Fig. ||). For sufficiently large values of n p , Auj <§; 7d 
even for 1^ = 1. 

4. UPPER ENVELOPE OF PARENT MODE AMPLITUDES 

Parametric instability provides an upper envelope to the 
amplitudes of overstable modes. Coupling of an overstable 
parent mode to a single pair of daughter modes suffices to 
maintain the pare nt mode's amplitude near the thresh- 
old value (cf. §2.2). The energy gained by the overstable 
parent mode and transferred to the daughter modes may 
be disposed of by linear radiative damping or by further 
nonlinear coupling to granddaughter modes. 

The results presented in this section are obtained from 
numerical computations and displayed in a series of fig- 
ures. The general trends they exhibit are best understood 
in terms of analytic scaling relations. We derive these first 
in order to be able to refer to them as we describe each 
figure. 

Excited modes of ZZ Ceti stars are usually detected 
through photometric measurements of flux variations, and 
in a few cases through spectroscopic measurements of hor- 
izontal velocity variations. Hence we calculate surface am- 
plitudes of fractional flux, AF/F, and horizontal velocity, 
Vh, variations. 9 Each of these is directly related to the 
near surface amplitude of the Lagrangian pressure pertur- 
bation, Sp/p. The threshold value of the latter is obtained 
by combining equations (Q) and (12) with the normaliza- 
tion factor (nT^L)^ 1 / 2 given by equation (A28) of Paper 
I. Thus 



Sp 

p 



Adopting relations expressing vu and AF/F in terms of 
Sp/p from §3 of Paper I, we arrive at 



LuRn 



Vh 



1/2 



(17) 



and 



AF 
~F~ 



1 + K)' 



1/2 



(1±\ + ( °J± 

\ U) J \ UJ 



SuY 
uj J 



1/2 



(18) 



In the above, R denotes the stellar radius, and r c is the 
thermal time constant describing the low pass filtering ac- 
tion of the convection zone on flux variations input at its 
base. 10 

Figure ^| displays calculated values for amplitudes of 
overstable modes limited by parametric instability. The 

9 In this section we drop subscripts on parent mode parameters and 
10 t c ss 3-rth, where Ty, is evaluated at Zf,. 



rise of with increasing mode period mainly reflects 

the correspondi ng r ise of the damping rates of the daugh- 
ter modes (cf. §|3.2| ) as indicated by equation (|l6|). At the 
longer periods, the values of |<5p/p| decline with decreasing 
T e ff. This is a subtle consequence of the deepening of the 
convection zone which pushes down the top of the daugh- 
ter modes' cavities, thus redu cing their damping rates (see 
§5.2 of Wu & Goldreich |l999| , hereafter Paper II). The be- 
havior of \vh\ is similar to that of |<5p/p| except that \vh\ 
decreases relative to |<5p/p| with increasing £ as shown by 
comparison of equations ([l6]) and jl7). A new feature 
present in the run of |<5F/Fjverses mode period is the low 
pass filtering action of the convection zone as expressed 
by the factor [1 + (wt c ) 2 ] -1 / 2 in equation (18). This fac- 
tor causes \SF/F\ to rise slightly more steeply than |<5p/p| 
with increasing mode period. It is also responsible for a 
more dramatic decrease in \8F/F\ with decreasing T e g at 
fixed mode period. Since \vh\ does not suffer from this 
visibility reduction, velocity variations may be observable 
in stars that are cooler than those at the red edge of the 
instability strip. 

Figure @ reproduces a summar y of o bservational data on 
mode amplitudes from Clemens (1995). Each star is repre- 
sented by a point showing the relation between the V-band 
photometric amplitude in its largest mode and the ampli- 
tude weighted mean period of all its observed modes. The 
latter quantity is a surrogate for the star's effective tem- 
perature in the sense that longer mean per iods c orrespond 
to lower effective temperatures (Clemens 1995, Paper I). 
Our theoretical predictions for the amplitudes of £ = 1 
modes limited by parametric instability are shown by a 
solid line obtained by interpolation from the numerical re- 
sults displayed in Figure [5]. Two features of this figure are 
worthy of comment. 

For a few low order overstable modes, n p < 3 at £ p = 1, 
Suj is more significant than 7^ in determining the best 
daughter pair (Fig. ||). The best daughter pairs for these 
modes have £ values of a few and identities which depend 
sensitively on minor differences among stars. Thus we ex- 
pect amplitudes of short period overstable modes to show 
large star to star variations. That the n = 1,£ = 1 mode 
is detected in only about one half of the hot DAVs (see 
Figure 4 of Clemens 1995) is consistent with the expected 
statistical variations in Suj. 

Overstable £ = 1 modes with periods longer than 800 s 
(n p > 20) are probably not saturated by parametric in- 
stability. Their maximum |k| is severely reduced below 
the adiabatic value by t he strong nonadiabaticity of their 
daughter modes (§B.2). Other mechanisms that con- 
tribute towards saturating these modes are discussed in 



5. DISCUSSION 

5.1. Turbulent Saturation 

Turbulent convection severely reduces the vertical gradi- 
ent of the horizontal velocity of g-modes in the convection 
zone. As a result, a shear layer forms at the boundary be- 
tween the bottom of the convection zone a nd the top of the 
radiative interior (Goldreich & Wu 199£, hereafter Paper 
III). Kelvin-Helmholtz instability of this layer provides a 

denote daughter mode parameters by a subscript d. 
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Fig. 6. — Selection of the best daughter pairs for a few low order £ p = 1 parent modes in a white dwarf model with T c g = 12, 800 K. We 
look at daughter pairs with £ dl ranging from 1 to 100, and l d = £ d + £ v . We plot (7d 1 7d 2 ) 2 m the upper panel, and in the middle 
panel, for the pair that sets the lowest instability threshold at each £ dl . This does not necessarily correspond to the pair that give rise to the 
lowest \Su>.\ at each £ dl . The corresponding threshold value of (Sp/p) p evaluated at the stellar surface is shown in the bottom panel. With 
increasing £ dl , ("/d 1 ld 2 ) 1/l2 rises and \6ui\ becomes unimportant. The best daughter pairs have low £ values. 



nonlinear dissipation mechanism for overstable modes. An 
overstable mode's amplitude cannot grow beyond the value 
at which nonlinear damping due to the Kelvin-Helmholtz 
instability balances its linear convective driving. Equa- 
tion (44) of Paper III provides an estimate for the value 
at which this mechanism saturates the surface amplitude 
of (Sp/p); 



0.1 [(^r c ) 2 + l] 1 /2[( WTc )2_ 1]L 



C D 



0JT c 



Rzb 



(19) 



Our ignorance of the complicated physics involved in a 
nonlinear shear layer is covered by the range of possi- 
ble values of the dimcnsionless drag coefficient, Gr>. Ter- 
restrial experiments indicate that Co falls between 10 -3 
and 10 _1 . The dashed lines in Figure || show the ef- 
fect of including nonlinear turbulent damping in addi- 
tion to parametric instability on limiting mode amplitudes. 
Amplitudes of overstable modes saturated by the Kelvin- 
Helmholtz instability are uncertain because Cd is poorly 
constrained. 



be close to their equilibrium values. This raises a worry. 
Suppose the daughters are prevented from reaching their 
equilibrium amplitudes by parametric instability of grand- 
daughter modes Then they would not be able to halt the 
amplitude growth of the parent mode. 

To answer the first of these questions, we calculate the 
ratio, denoted by the symbol S, between the threshold 
amplitude for a daughter mode to excite granddaughter 
modes and its equilibrium amplitude under parametric ex- 
citation by the parent mode. The former is obtained from 
equation (Q), and the latter from equations (l|)-(R). We 
make a few simplifying assumptions to streamline the dis- 
cussion. Resonances between daughters and granddaugh- 
ters are taken as exact; individual members of daughter 
and granddaughter pairs are treated as equivalent. Equa- 
tions (|l2|) and (Ful) are combined to yield 



,,2 ■ > , 

n 2 p L 

It is then straightforward to show that 



(20) 



5.2. Granddaughter Modes 

Here we answer the following questions. Under what 
conditions do daughter modes excite granddaughter modes 
by parametric inst ability? 11 What are the consequences if 
they do? 



In 



3 2.2 we show how parametric instability of 
linearly damped daughter modes maintains the amplitude 
of a parent mode close to its equilibrium value. In order to 
dispose of the energy they receive from the parent mode, 
the time averaged energies of the daughter modes must 

In this subsection subscripts p, d, and g refer to parent, daughter. 



LUpUj d / n d 



7, 



n p J \ 7d + 



32 



12 + ^ 2 



(21) 



The factor 32 is an approximation based on taking 
ujg/ujd = Ud/up = 1/2 and n p /n d = 1/2. 

In general S ^> 1, so the excitation of granddaughter 
modes requires the daughter modes to have energies in ex- 
cess of their equilibrium values. However, the equilibrium 
solution is unstable if the best daughter pair corresponds 

and granddaughter modes. 
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Fig. 7. — An assembly of results on threshold amplitudes for parametric instability. From left to right, amplitudes of overstable modes 
with different £ are plotted against mode period. From top to bottom, the panels display photospheric amplitudes of fractional pressure 
perturbation, horizontal velocity, and fractional flux variation. Filled circles, crosses, and open squares symbolize modes from the three stellar 
models in Figure a in order of descending T c ff . Pulsation amplitudes generally rise with increasing mode period, except for the shortest 
period modes which exhibit large star to star fluctuations. Amplitudes of \Sp/p\ for I = 1 modes in the coolest model dip around 600 s and 
flatten beyond 1200 s. The first feature is associated with the dip in 7 for neutrally stable daughter modes, and the second is explained by 
the flattening of 7 for strongly nonadiabatic daughter_xnodes of long period (see Fig. H). Reduction of the coupling coefficient due to strong 
nonadiabaticity is not accounted for (see appendix p.2|) . 



to 7^ > \Su\, and then the daughter mode energies episodi- 
cally rise far above their equilibrium values. At such times, 
granddaughter modes may be excited by parametric in- 
stability and consequentially limit the amplitude growth 
of the daughter modes. This slows the transfer of energy 
from parent to daughter modes, but it does not prevent 
the daughter modes from saturating the growth of the par- 
ent mode's amplitude at the level described by equation 

For the few lowest order parent modes, we typically find 
1 5a; I > 7<i. This may reduce S to below unity with the con- 
sequence that the parent mode amplitude may rise above 
that given by equation (^). 

5.3. Additional 3- Mode Interactions 

Parametric instability sets reasonable upper bounds on 
the photospheric amplitudes of overstable modes. In a 
given star this upper bound rises with increasing mode pe- 
riod except possibly for the lowest few modes. However, 
the observed amplitude distributions are highly irregular. 
This mode selectivity may arise from 3-mode interactions 
which involve more than one overstable mode. 

We investigate a particular example of this type. It 
is closely related to parametric instability, the only dif- 
ference being that the daughter modes of the overstable 
parent mode are themselves overstable. Acting in isola- 
tion, resonant mode couplings tend to drive mode ener- 
gies toward equipartition. They conserve the total energy 



E p + Ed t + E c i 2 — 0, and transfer action according to 
E^/wd-L — Ed 2 /oJd 2 = —E p /lu p . In this context it is im- 
portant to note that the energies of overstable modes lim- 
ited by parametric instability decline with increasing mode 
period (Fig. ph. Therefore nonlinear interactions transfer 
energy from the parent mode to its independently excited 
daughters. As shown below, this transfer may severely 
suppress the parent mode's amplitude. 

We start from equations ([!])- (|3|). These may be manip- 
ulated to yield 

dE 

—L = lp E p + 3V2LOp K {EpE dl Ed 2 ) 1/2 sin$, (22) 

where $ = 9d 1 +0 d2 —6 p . For E p 3> E dl and Ed 2 , nonlinear 
interactions transfer energy from the parent mode to its 
daughter modes. In particular, if we ignore phase changes 
in the overstable daughter modes due to their interactions 
with granddaughter modes, we find that $ satisfies 



E, 



, 1 /2 



E dl 



Ud 2 
E d2 



E„ 



cos $, 
(23) 



with a stable solution at $ = — 7r/2 when 8lj = 0. 

We denote the ratio of the nonlinear term to the linear 
term in equation (p3) by the symbol T; 



T 



3V2a 



7 P 



Ed x Ed 
E n 



1/2 



(24) 
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Fig. 8. — Comparison between theoretical parametric amplitudes and those summarized by Clemens (1993) from observations of two dozen 
DA variables. Each star is represented by a filled circle plotted such that the ordinate gives the photometric amplitude of the largest mode 
and the abscissa the mean period weighted by the photometric amplitudes of all its detected modes. To compare theory and observation, we 
single out the largest amplitude I = 1 mode for each of the three models used in Figure M. The solid line is an interpolation between these 
largest amplitude modes and their periods. The amplitudes are reduced by a factor of 2.5 to approximately account for the reduction in visible 
relative to bolometric flux variations for DA variables. The run of observed mode amplitude as a function of mode period is mimicked by the 
solid line. However, the extrapolation to short period shown by the dotted line is not reliable, because in this regime mode amplitudes have 
large star-to-star variations associated with different values of frequency mismatches. The drop-off in amplitude at long period (P > 10 3 s) 
results from a combination of lower daughter mode damping rates and a reduction in the visibility of the parent mode due to increasing u>t c . 
Turbulent dissipation arising from the shear layer below the convection zone (§ p,l| ) further suppresses the amplitudes that long period modes 
can attain. This is illustrated by the dashed curves for different choices of the drag coefficient, Crj>. 



Using the magnitudes of Ei set by parametric instability 
of their respective dau ghte rs, and adopting the same ap- 



proximations made in §5.2 



T 



Td + 5uj2 



we arrive at 



32 



Id + 5uj2 



(25) 



Comparing equations ( |2l| ) and ( |25| ) we see that S = T . A 
little thought reveals that this is not a coincidence. For 
T>1, overstable daughter modes can suppress a parent 
mode's energy below the value set by parametric instabil- 
ity. We expect this suppression to be important in cool ZZ 
Ceti stars whose overstable modes extend to long periods. 
It may render their intermediate period modes invisible. In 
a similar manner, the amplitudes of high frequency over- 
stable modes with i — 2 and 3 may be heavily suppressed 
by interactions with their lowest I overstable daughters. 

The irregular amplitude distributions among neighbor- 
ing modes may be partially accounted for by this type of 
resonance. Mode variability may also play a role. We 
explore this in the next subsection. 

5.4. Mode Variability 

Excited g-modes in ZZ Ceti stars exhibit substantial 
temporal variations. Parametric instability may at least 
partially account for these variations. 

When \Suj\ < parametric instability gives rise to 
limit cycles in which the amplitudes and phases of parent 



and daughter modes vary on time scales as short as 7J 1 - 
Stable daughter modes may briefly attain visible ampli- 
tudes. Temporal amplitude variations may contribute to 
the irregular mode amplitude distribution seen in individ- 
ual stars. 

Phase variations of a parent mode obey the equation 



d6 p 
~dT 



3 \A dl \\Ad 2 



cos $. 



(26) 



At the equilibrium given by equations (H)-®, t ne parent 
mode's frequency is displaced from its unperturbed value 
such that 



d9 p 
~dt 




This constant frequency shift is of order 10~ 9 s~ 1 for the 
n = 1, £ = 1 mode and of order 10 -7 s -1 for the n = 2, 
1 = 2 mode. Frequency shifts in higher order overstable 
modes which are involved in limit cycles are predicted to be 
larger and time variable. During brief intervals of length 
~ 7rf' wnen the daughter mode energies are comparable 
to that of the parent mode, \to p — w p \ ~ 7<j, which is of order 
a few times 10 _5 s _1 , or a few 11Hz in angular frequency. 
These shifts might account for the tim e- vary ing rotational 
splittings reported by Kleinman et al. (199?) provided dif- 
ferent m components of the overstable modes are involved 
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Fig. 9. — Energy versus period for overstable modes in the coolest model considered in Figure W. Solid, dashed and dotted lines represent 
I = 1, 2 and 3 modes. Mode energy tends to decrease with increasing period and with increasing t. Deviations from this trend for the I = 1 
modes are explained as follows. The narrow dip near 200 s is due to a fortuitiously good frequency resonance involving the n = 2 mode, and 
the wider dip around 700 s is associated with neutrally stable daughter modes. 



in different limit cycles. 

5.5. Miscellany 

We briefly comment on two relevant issues. 

Gravitational settling produces chemically pure layers 
between which modes can be partially trapped. Modes 
that are trapped in the hydrogen layer have lower mode 
masses, and therefore higher growth rates and larger max- 
imum values of \k\ than untrapped modes of similar fre- 
quency. This implies lower threshold energies for para- 
metric instability. Nevertheless, trapping does not affect 
predicted photospheric amplitudes of Sp/p, and hence Vh 



and 5F/F, since these are proportional to A divided by 
the square root of the mode mass. 

In circumstances of small rotational splitting, the simple 
limit-cycles depicted in Figure [I] are unlikely to be realis- 
tic. In such cases, different m components of an overstable 
parent mode share some common daughter modes. This 
leads to more complex dynamics. 

We are indebted to Bradley for supplying us with mod- 
els of DA white dwarfs. Financial support for this research 
was provided by NSF grant 94-14232. 



APPENDIX 



THREE-MODE COUPLING COEFFICIENT 

We reproduce the expression for the three-mode coupling coefficient presented in equation (^|) of §3.1 with one modifi- 
cation: 

r ' (r ' 2) -(V-£) 3 + ~r 1 (v-oe 7 4 + Detie: 



where 



d 3 xp 
Det|^ 



G 



6^ ^ 2 



(V • + kxu* 



(Al) 
(A2) 



Because of strong cancellations among its largest terms, this form is not well-suited for numerical evaluation. We derive 
a new expression which does not suffer from this defect. Then we estimate the size of k and deduce its dependences upon 
the properties of the three-modes. 

Simplification 

In spherical coordinates (r, 9, </>), the components of the displacement vector may be written as 

&(r) 8 &(r) 1 d 



r dO 1 r sin 2 6 dcf) 



(A3) 



where Yi m is a spherical harmonic function, and the £j are covariant basis vectors. 
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The angular integrations in equation (Al) are done analytically. The following definitions and properties prove useful: 



T : 



dO sm6Y a Y b Y c =< Y a Y b Y c >, 



T 



F a = < Y a V a Y b V a Y c >= (A 2 + A 2 - A 2 a ) - 

G a = < v Q v >y v°nv /, y c >= [A* - (A 2 b - a 2 ) 2 ] |, 

S = < V a (V a Y a V p Y b V p Y c ) >= G a + G b + G c 



1 



-(A 2 ^ + A b F b + A 2 C F C ) = G a + A 2 a F a , 



2 

v a = < y a v Q v p nv /3 v Q y >= a 2 a 2 t -f*-s, 



A 2 a T = F b + F c , 
A 2 F a = G b + G c . 



(A4a) 
(A4b) 
(A4c) 

(A4d) 

(A4e) 
(A4f) 
(A4g) 



Here, V a is the covariant derivative on a spherical surface; a can be either 9 or cf>. Each angular integration is 
proportional to T which contains all m dependences and is of order unity independent of the i values of the participating 
modes. The paramter A 2 = £(£ + 1). Subscripts a, b and c denote different modes. The angular selection rules are simply 
4 e [14-41,4+4], Mod[(4+4 + 4),2] = and m a + m b + m c = 0. 

Angular integration of Det|£?-| yields 



tfflDet|&| = \ 



T£C - AlT&ft + \(A 2 A 2 C T - V a )&$ 



dtf 
dr 



Symmetrizing this expression with respect to modes b and c, employing equations ( A4d)-( A4g| ), we get 



(F a + S) d 

2r 2 dr 



A 2 T d T d 



2r 2 dr 



For gravity-modes, |^| 3> |£ r |. Thus the four terms in this expression decrease in size from left to right. 
In a similar manner we arrive at 



jdn (v • = nv ■ r)(v • £ 6 )(v • r), 



(A5) 



(A6) 



(A7) 



and 



—(v • £°)-r $8) + — (v • + nv • 

r dr r dr dr dr 



2TA 2 ,_ t „. . 2T. 



-(v-nea + ^i-(v-o^e; 



(A8) 



The magnitudes of the terms in the above expression decrease from left to right except that terms 2-4 are of comparable 
size. 

Having disposed of the angular dependences in n, we turn to the radial integrations. The following relations prove 
helpful in this context: 



dr r r 

d rr itt cm A V? t f 2 , 2 9 d 9 
— Tip(V ■ £) = 6i - w H — 

dr r \ r dr 

Tip(V ■ £) = pg£, r - rpu 2 £ h . 



(A9) 

(A10) 
(All) 



The first equation is the definition of divergence, and the second and third are the radial and horizontal components of 
the equation of motion written in terms of the Lagrangian displacement. 

Curvature terms arise because the directions of the basis vectors depend upon position. For instance, 



dY 



I III 



de ' 



(A12) 
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where £/, appears due to curvature in the coordinate system. The largest terms in equations ( |Aq) and (A8) are curvature 
terms. H owe ver, their radial integrals cancel leaving a much smaller net contribution. Direct numerical evaluation of 
equation ( |Al| ) leads to unreliable results, so it is i mpor tant t o ca rry out this cancellation analytically. Accordingly, we 
integrate each term by parts and apply equations (A4d) and (|A9|) to obtain 



dr 



SP d j c a c b c c\ 



dr 



F a T lP 



(A13) 



This step leads to 



. ^ r ^ r 

dr dr 



ar 2 



{S-F a )pg 



C a C° C c 



2 dr 



(A14) 



Next we systematically eliminate radial derivatives of the displaceme nt v ector from the expression for k. W e inte grate 
by parts to disp ose o f d^h/dr and substitute for d£ r /dr using equation (A9). With the aid of equations (A4d)-(A4i), and 
using equation (All) to make the expression symmetric with respect to indexes b and c, we arrive at our final working 
expression for k, 



dr i - 



dg 



TVryri + l)p f „ &a ^„ ttwT7 ^ , u 2 a G a rp b 



-(v-n(v-o(v-r 



c a c° £ c 



A'Trl^™ Q 



££(v ■ e°)(v • r) - r(3r a + i) P (v • OtC 



+2TrT lP C(V ■ £ h )(V • C) + ^CtrC + \ [H - 'i^l - 3^ c 2 ) F a - (2u? h F h + 2co 2 c F c )] rp^h 



1 

+ 2 



(A15) 



All permutations of the 3 modes are to be included when evaluating this expression. However, the largest contribution 
from each term comes when b and c are radially similar daughter modes. For high order gravity-modes, the firs t five 
terms are of comparable size and much larger than the remaining four. Numerical evaluation of k using equation ( A15| ) 
does not suffe r from the errors arising from large cancellations and numerical differentiation that plague attempts using 
equation (Al). 

Order- of-Magnitude 

Here we estimate the maximum value that k can attain for parametric resonances involving a given parent mode. In 
so doing, we apply results derived in Papers I and II. These include: (1) the scaling relations 



IV- €1 



A% Cr 



kh£,h 



1 



in the evanescent region z < z u , and 



IV -£| 



i A2 



z w \ 2 A 

K z 



{n TuJ L) 1/2 '' 



(A16) 



(A17) 



in the propagating cavity z > z^,; (2) the fact that regions between consecutive radial nodes contribute equally to the 
following normalization integral 



\j drr*p(k*ti +&)=!. 



(A18) 



As ^ <§C £,h f° r g-modes, we have 



u 2 R 2 



, a 2,2 n' 1 r , , 1 ( z 

dz pA E, h ~ — ~ — / dz k z ~ — — 



(A19) 
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Here n' is th e num ber of radial nodes above depth z and n the total number in the mode. 

Equation (A15) yields maximal values for k when mode a is taken to be the parent mode and modes b and c to be 
two radially similar daughter modes. Most of the contribution to the radial integral comes from the region above z lAJa : 
for parametric resonance, z UJa is much greater than z Ub ~ z UJc ; the decay and rapid oscillation of the parent mod e's 
eigenfunction renders insignificant contribution from greater depths. Thus we can take the integrals in equation ( A15| ) to 
run from z — to z — z Ua and pull out the parent mode eigenfunctions since they are approximately constant for z < z tJa . 
This procedure reduces each of the leading terms in k and thus their sum to 



1/2 



ridinpT^L) 1 / 2 ' 



But z UJp /z 0Jd ~ A-l/A^ ~ n 2 d /n 2 p which leads to 



1 



(A20) 



(A21) 



The above equation establishes that the maximum value of k rises steeply with increasing radial order of the parent mode 
and is independent of the radial orders and spherical degrees of the daughter modes. 



STRONGLY NONADIABATIC DAUGHTER MODES 
Strong nonadiabaticity occurs wherever 

SSI- 



(k z z) 2 



(Bl) 



We refer to the depth above which this inequality applies as z na . To derive a simple analytic scaling relation for z na , 
we make use of the approximations k z ~ (zz w ) -1 / 2 and Tth/Tb ~ (z/z b ) q+2 . Here z& and Tb are the depth and thermal 
relaxation time at the bottom of the convection zone, and p oc z q with q « 3.5 provides a fit to the density structure in 
the upper portion of the radiative interior. It then follows that 



Zna 
Zb 



1 z b 

U)T b Zu 



Moreover, 



(Mr 



By reducing the effective buoyancy, strong nonadiabaticity lowers the effective lid of a g-mode's cavity to z na 
of this fact are explored in the following subsections. 

Damping Rates of Strongly Nonadiabatic Modes 
As shown in Paper II, the energy dissipation rate for a strongly nonadiabatic mode may be written as 

to 



(B2) 
(B3) 

Consequences 



7 



lnTT 



where TZ denotes the coefficient of amplitude reflection at the top of the mode's cavity. 

To derive an approximate relation for TZ, we note that the real and imaginary parts of k z , k zr and k z 



satisfy 



(k zr z) 2 



provided |& zi |/|fc zr | < 1. Thus 



lnTT 



dz \k z 



Evaluating this integral with the aid of equations (B2) and (B3), we obtain 

^ £(«+2)/(g+l) 



lnTT 



(Wfe) 



q+2 



w (2q+5)/(2g+2) 



(B4) 



(B5) 



(B6) 



(B7) 



Numerical results for In TZ 1 plotted in the upper panel of Figure BIO confirm this relation 



Because the maximum value of k for parametric resonance is independent of t he sp heri cal degrees of the daughter 
modes, the dependence of 7 on I at fixed u> is of great significance. Equations (B4) and (B7), together with the relation 
n oc I at fixed oj, establish that 7 oc Since 7 increases with I at fixed lo, the most important daughter pairs are 

those with the smallest £ values subject to the constraint 7^ > 6u>. 
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Fig. BIO. — Strongly nonadiabatic modes in the hottest model considered in Figure |5| The top figure exhibits the dependence of the 
reflection coefficient, 1Z, on £ for modes with periods close to 550 s. The numerical result for the reflection cocfficieni_is obtained from 
numerical computations of 7 and is depicted by filled triangles. The underlying solid line is the theoretical estimate (eq. EST]]). Computing 7 
requires care since 1Z s=s 10 -16 for i = 50. This is of order the machine accuracy for double precision numbers. The bottom panel shows the 
angular degree at which decoupling occurs, £<j c , as a function, of period for 1 = 1 parent modes. Triangles denote numerical values and the 
solid line the analytical scaling relation given by equation (pa). We find that i& c is almost independent of the star's effective temperature. 



Reduction of the Coupling Coefficient by Strong Nonadiabaticity 

The maximum adiabatic coupling coefficient between a parent mode and a pair of daughter modes is K m ax ~ 
(^pT Wp i) -1 / 2 (cq. fl^| ). Its major contribution comes from the region above z UJp . The factor (n p T Up L)~ 1 / 2 is the 

surface value of the normalized eigenfunction for [£(( + l)] 1 / 2 ^/ R. The extra factor n" 1 is the fraction of each daughter 
mode's nodes that lie above z UJp . The coupling coefficient is reduced compared to equation ( |l2|) for daughter modes that 
arc strongly nonadiabatic in the region above z u]p . 

Equation (Bl) indicates that nonadiabaticity increases with increasing I at fixed z and uj; k z ~ (zz w ) -1 / 2 and z w ~ 
uj 2 R 2 I ' gl{l + 1), so k z oc £. We define the angular degree of decoupling for a given parent mode, 4o as the smallest 
spherical degree at which its daughter modes are s tro ngly nonadiabatic all the way down to the top of the parent mode's 
cavity; that is, z nad w z u at i^c- Using equation (B2), we find 



4c 



gzb 



LUT b 



1/2 



(B8) 



Numerical results displayed in the lower panel of Figure B10| are well-represented by this scaling relation. We find that 
£dc is relatively independent of stellar effective temperature but decays steeply with mode period. By n p = 20, we find 

4c < 1. 

It is plausible that at id = 4c, K max is reduced by a factor n p /nd below its adiabatic value (eq. pH ) because the 
effective lids of the daughter modes' cavities are lowered to z UJp . Since n p /na ~ £ p /2lg, this is a large reduction for 
4c 3* 4- A n even more severe reduction of K max is expected for t& > Iso- 
parametric Instability for Traveling Waves 

In the limit of strong nonadiabaticity, daughter modes are more appropriately described as traveling waves than as 
standing waves. Thus it behooves us to investigate the parametric instability of traveling waves. Here we demonstrate 
that the instability criterion for traveling waves is equivalent to that for standing waves (eq. |Q). 

Nonlinear interactions between parent and daughter modes are localized within an interaction region above z UJp . Let 
us assume that z nad <C z UJp . Then, in most of the interaction region daughter wave packets may be represented as linear 
superpositions of adiabatic modes. Propagating at their group velocity, the daughter wave packets pass through the 
interaction region in a time interval 

dz 1 Tm r i 

-. (B9) 



AT = / 
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Three significant relations involving n p are worth noting: 1) n" 1 is the fraction of each daughter mode's nodes that 
lie above z UJp , so 2 AT is a fraction n~ 1 of the time each daughter wave packet takes to make a round trip across its 
cavity; 2) approximately n p daughter modes reside within the frequency interval n/ AT, and their relative phases change 
by less than 7r as each daughter wave packet crosses the interaction region; 3) maximal k occurs inside an interval of width 

|rc«2i ~ n d 2 \ £ n V 

Nonlinear interactions between parent and daughter waves within the interaction region are described by equations 
(0)"(H) w ith two modifications of the equations governing the time evolution of the daughter modes. The linear damping 
term is negligible for z 3> z nai and the nonlinear term must be multiplied by a factor n p . The latter accounts for the 
number of modes which couple coherently to each daughter mode during the interaction time AT. During two passes 
through the interaction region, the amplitudes of the daughter wave packets grow by a factor e G , where the gain, Q, is 
given by 

G= 2 ^^*3V2n d \K\\A P \. (BIO) 

For parametric instability to occur, 

G>lnTZ- 1 . (Bll) 



Combining the relation between 1Z and 7 given by equation (B4) with equation (Bll), the threshold condition for 
parametric instability of traveling waves becomes 

' P ' 3V2u> d \K\ 

The above condition is equivalent to equation (|J) in the limit that 7^ 3> Su>. Thus the threshold condition for parametric 
instability of traveling waves reduces to a limiting case of the threshold condition for parametric instability of standing 
waves. 
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